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ABSTRACT 

In high-redshift {z > 5) quasars, a large amount of dust (r^ 10^ Mq) has been observed. 
In order to explain the large dust content, we focus on a possibility that grain growth by 
the accretion of heavy elements is the dominant dust source. We adopt a chemical evolution 
model applicable to nearby galaxies but utilize parameters adequate to high-z quasars. It is 
assumed that metals and dust are predominantly ejected by Type II supemovae (SNe). We 
have found that grain growth strongly depends on the grain size distribution. If we simply use 
the size distribution of grains ejected from SNe, grain growth is inefficient because of the lack 
of small grains (i.e. small surface-to-volume ratio of the dust grains). However, if we take 
small grain production by interstellar shattering into consideration, grain growth is efficient 
enough to account for the rich dust abundance in high-z quasars. Our results not only confirm 
that grain growth is necessary to explain the large amount of dust in high-z quasars, but also 
demonstrate that grain size distributions have a critical impact on grain growth. 

Key words: dust, extinction - galaxies: evolution - galaxies: high-redshift - galaxies: ISM - 
ISM: clouds - quasars: general 



1 INTRODUCTION 



The origin of dust grains in the Universe is one of the funda- 
mental problems in astrophysics. It is widely believed that the 
evolution of dust content in galaxies is influenced by dust for- 
mation in stellar ejecta, dust destruction in supernova (SN) rem- 
nants, and grain growth by the accretion of metals onto t he preex- 
isting dust grains in molecular clouds (e.g. lDwek|[l998b . Various 
authors have shown that grain growth in molecul ar clouds domi- 
nates t he inc r ease of dust mass in nearby galaxies dHirashitall 19991 : 
Inoue 2003; Zhukovska et al. 2008!; iDrainel |2009| ; llnoud 1201 ll) . 
Zhukovska et aL ( 2008) also mention that grain growth depends on 
the mean surface-to-volume ratio of the dust grains. In our previ- 
ous work l lHirashita&Kuoll201ll. hereafter HKl 1), we have inves- 
tigated the effects of dust grain size distribution on grain growth, 
and found that the grain size distribution has a large influence on 
the efficiency of grain growth because the surface-to-volume ratio 
is governed by the grain size distribution. 

The extinction curves in the Galaxy and the Magellanic 
Clouds can be explained with dust grain size distributions n(a) oc 
a~^'^, where a is the grain radius and n{a)da is the num- 
ber d ensity of grains whose radius is between a and a + da 
jMathis. Rumpl. & Nordsieck 1977; Pei 1992). However, high- 
redshift quasars ex hibit different extinctio n curves from the ones in 
local galaxies ("e.g. iGallerani et al ] l2O10h . which may indicate dif- 
ferent grain size distributions. Different dust sources and/or differ- 
ent grain processing in the early Universe may be the reason for the 



different extinction curves in high-z quasars. In the early Universe, 
Type II supemovae (simply denoted as SNe in this paper) should 
be the foremost source o f dust formation because of the short life- 
times of the progenitors et al.lbOlll) . albeit asymptotic giant 
branch (AGB) stars may ha ve contribution even at z > 5 for some 
stellar initial mass function jValiante et al.l l2009'). Star-forming re- 
gions may also host s ignificant grain processing such as shattering 
( iHirashita eTai]|2010l) . 

Grain growth in molecular clouds may dominate the grain 
mass i ncrease not on ly in the local Univers e but a lso at 
high z jMattssonl 1201 ll : IPipino et all 1201 ll ; IVahante et alj I2OI ll ; 
lAsano et al .I2OI2I) . However, the potential importance of dust grain 
size distributions in grain growth has not been focused on. As 
emphasized above, the grain size distributions in high-z quasars 
should be different from those in local galaxies. Therefore, we 
should reconsider grain growth by using grain size distributions ap- 
plicable to high-2 systems. 

Size distributions of grains produced and ejected from SNe 
(called SN-dust) are calculated bv lNozawa et al.l ( l2007h . who take 
dust condensation and destruction in SNe into account. They show 
that the average grain size of SN-dust is biased to large sizes be- 
cause small grains are efficiently trapped in the shocked region and 
considerably destroyed by thermal sputtering. This means that the 
grain size distribution is completely different from the local a~'^'^- 
law and that the surface-to-volume ratio is smaller than the local 
one. As a result, we have to reexamine whether grain growth is effi- 
cient enough to explain the large dust abundance in high-2 quasars. 
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After SN-dust is injected into tiie interstellar medium (ISM), 
tlie grain size distribution can be modified. Dust grains acquire 
large velocity dispersions in a warm ionized mediu m through dy- 
namic al coupling with turbulent motion ( Yan. Lazarian. & Draind 
l2004h . Grain-grain collisions at velocities driven by turbulence 
lead to efficient shattering, w hich produces a large number of small 
grains jHirashita et al. 2010h. Th e extin ction curve observed in a 
high-2 quasar bv ' Maiolino et alj (|2004) is also reproduced by a 
grain size distribution after shattering (Hirashita et al . 2010) . Shat- 
tering raises the surface-to-volume ratio, leading to a higher grain 
growth rate. 

In this Letter, we examine the significance of grain size dis- 
tribution on grain growth in high-z quasars by focusing on two 
different grain size distributions as representative ones expected 
theoreticall y in high-z quasars: one is the grain size distribution 
of SN-dust jNozawa et all2007l) , and the other is the size distribu- 
tion of grains processed by interstellar shattering dHirashita et al.l 
|20l3). This Letter is organized as follows. In Section 2, we explain 
theoretical models of dust enrichment, focusing on grain growth 
by accretion. In Section 3, we show the results, which we discuss 
in more general contexts in Section 4. Finally, Section 5 gives the 
conclusion. 



2 MODEL 

2.1 Dust Enrichment 

We adopt a chemical evolution model that describes the time evo- 
lution of gas, metals, and dust in a galaxy. In order to focus on the 
effect of grain size distribution on grain growth, we adopt a simple 
model, which has alread y been shown to be applicable to nearby 
galaxies l lHirashitall999l : HKll), but use parameter values suitable 
for high-z (z > 5) quasars. We neglect inflow and outflow, and 
treat the system as a closed boxQThe abundances of dust and met- 
als are tightly related, so the model equations are finally reduced to 
the relation between dust-to-gas ratio and metallicity by adopting 
the instantaneous recycling approximation (HKll), which can be 
applied as long as we consider time-scales longer than the lifetime 
of massive stars: 



/i„,x(7^Zx + Vx) - WsN + 7^)Dx 
dMdust.x" 



1 



dt 



(1) 



where Zx is the mass abundance of element X (X is a key element 
of a dust species; X = Si for silicate and X = C for carbonaceous 
grains) in both gas and dust phases, Ox is the mass abundance 
of a dust species whose key element is X, TZ is the returned frac- 
tion of the mass from stars formed, 3^x is the mass fraction of el- 
ement X that is newly produced and ejected by stars, /in,x is the 
dust condensation efficiency of the element X in the ejecta, ^ is 
the star formation rate, j3sN is the destruction efficiency of dust by 
SN shocks (defined below), and [dA/dust,x/df]acc is the increas- 
ing rate of dust mass by accretion. The destruction efficiency by 



^ As long as the inflow has much smaller abundance of metals and dust 
than the system, both the metallicity and the dust-to-gas ratio decrease with 
the same fraction. Moreover, as long as the mass loss, for instance, by mass 
accretion onto the central black hole (BH) or outflow of gas heated by the 
feedback from BH or stars, has the same metallicity and dust-to-gas ratio 
as the host galaxy, it does not affect the relation between metallicity and 
dust-to-gas ratio. 



SN shock was introduced bv lMcKed h989l) as /3sn = e^M.-y/i;, 
where eg is the fraction of dust destroyed in a single SN blast, and 
A/s is the gas mass swept per SN blast, and 7 is the SN rate. We 
adopt PsN ~ 9.65 (HKll). In deriving equation (T]l, some quan- 
tities are eliminated. The gas mass is used to normalize the metal 
and dust masses to obtain Zx and Dx, respectively, and thus acts 
as a normalizing factor in our formulation. The star formation time- 
scale is also eliminated: changing the star formation time-scale just 
varies the quickness of metal and dust enrichment, leaving the re- 
lation between Zx and Dx unchanged. We refer to HKll for the 
derivation of equation ([TJ. We assume Dx = at Zx = 0. 

2.2 Grain Growtii by Accretion 

In HKl I, we have formulated the grain growth by accretion consid- 
ering the dependence on the grain size distribution. We also assume 
that the molecular clouds host both grain growth and star formation. 
Then, we obtain (HKll) 



dAfd 



dt 



(2) 



where j3x is the mass increment of dust whose key species is X and 
e is the star formation efficiency in a single molecular cloud. By 
assuming that grain growth is governed by the sticking of the key 
element X (i.e. Si for silicate and C for carbonaceous dust), HKI 1 
derived the following formula for fix'- 



Px 



Cx 



3y(a2)o + 3y2(a)o + j/3 



Cx 



(3) 



where y = ao^xTci/T (o-o = 0.1 [im is a typical grain radius, t^i 
is the lifetime of molecular clouds, and r oc ao is the grain growth 
time-scale given below in equation ID note that ao cancels out in 
obtaining y so that we can take an arbitrary value for ao), (a*)o 
is the average of (£ = 1, 2, and 3) weighted for the grain size 
distribution (i.e. moments), and fx = 1 — is the fraction 

of element X in gas phase (fx is the mass fraction of element X 
contained in the grains). The grain growth time-scale r is given by 



r = ao 



mxSx 
fxpx 



^x0 



(i)e"«( 



27rmx 



(4) 



where nu is the number density of hydrogen nuclei in the molec- 
ular clouds, (X/H)q is the solar abundance of element X relative 
to hydrogen in number, px is the material density of a dust grain 
whose key species is X, Sx is the sticking probability of element 
X onto the dust, mx is the atomic mass of element X, Tgas is the 
gas temperature in the molecular clouds, and is the Boltzmann 
constant. By using equation equation (TJ can be rewritten as 



=/in,x(7^Zx+3^x) 

dZx 



2.3 Grain size distribution 



iPsN + n- l3x/e)Vx. 



(5) 



We assume that the dust and metals are predominantly supplied 
from SNe in high-z (z > 5) quasars because they are the first 
significant dust source (Introduction). Although AGB stars may 
contribute to the dust enrichment even at 2 > 5, our conclusion 
in this paper is not altered as long as AGB star s supply large 
(~ . 1 ^.m) grains as su ggested by observations ( IGroenewegenl 
ll997l : lGaugeretal.lll999h . The dust and metal yields are repre- 
sented by the cas e of 20 M^d progenito r with an unmixed helium 
core, according to lHirashita et al.l(l2010l) (dust and metal yields are 



Dust growth in high-redshift quasars 3 



taken from iNozawa et alj I2OO7I and lUmeda & Nomotol I2OO2I re- 
spectively). The grain species formed by SNe are C, Si, Si02, Fe, 
FeS, AI2SO3, MgO, MgSiOa and Mg2Si04. Since the chemical 
properties of these species in grain growth is poo rly known, we 
simpl y categorize them into two groups following iHirashita et alJ 
( 120 IQl : the carbonaceous species and the silicate species (all the 
species other than C; Si is dominated), and hypothesize that they 
evol ve separately a nd independently in molecular clouds. 

iNozawa et all 12007.) have calculated the grain size distribu- 
tion by considering dust condensation and destruction in a SN, 
and have shown that small grains are efficiently trapped in the 
shocked region and considerably destroyed by thermal sputtering, 
and thus the size distribution of grains finally ejected into the 
ISM is biased towards larger sizes. Because destruction of small 
grains with a < a few x 10^"^ (i.m is common regardless of the 
progenitor mass and the treatment of mixing in the He core in 
[Nozawa et al. (2007) (similar grain size distributions are also ac- 
quired by IBianchi & Schneidedl2007h . it is sufficient to examine 
a single representative case for the grain size distribution for SN- 
dust. This case is ca ll ed 'gr ain size distribution without shattering'. 

IHirashita et aD ( I2OIOI) have suggested that grains are pro- 
cessed by shattering after being injected into the ISM. Production 
of small grains by shattering increases the total surface of grains, 
which is expected to raise the grain growth e fficiency. We adopt the 
grain size distribution after shattering from IHirashita et alj ^20 id) 
(the solar-metallicity case with hydrogen number density tih = 1 
cm~'^ and shattering duration t — h Myr). This case is called 'grain 
size distribution with shattering'. Although it is difficult to con- 
strain t and wh from observations, this single case is sufficient to 
demonstrate the importance of the change of grain size distribu- 
tion by shattering. Since we only focus on grain growth around 
Z > Zq, the assumption of solar metallicity is reasonable in this 
Letter. 

The above two grain size distributions adopted are shown in 
Fig.Ejfor silicate and carbonaceous species. The moments of grain 
radii {{a^)o) necessary to estimate /3x in equation ([3} are listed in 
Table [T] We observe that shattering makes the typical grain size 
smaller. We fix (a*)o; in other words, we fix the grain size distribu- 
tion throughout the galaxy evolution to focus on the effects of grain 
size distribution on grain growth. 

2.4 Selection of parameter values 

The typical time-scale of grain growth in equation (|4j is estimated 
as (HKIl): 

r = 6.30 X lO^ao.i (Zsi/Zg;^©) n3"'r5o"'/'5'o.3"' yr (6) 
for silicate, and 

r = 5.59 X lO'^ao.i (Zc/^c,©) 



-Irr -1/2 
J 50 



S0.3" yr (7) 

for carbonaceous species, where ao.i = ao/0.1 |J.m, ns = 
r^H/10^cm-^ Tgo = rgas/50K, and 50.3 = Sx/O.S. We 
adopt the same values as HKIl (see also the references therein): 



riH 



lO'^ cm-^ T^, 



50 K, and S'x = 0.3. For the solar abun- 



dance of Si and C, we adopt (X/H)q listed in Table|2l and convert 
it to Zx by Zx = (mx/^mH)(X/H) — 1.4 is the coiTec- 
tion for helium). We adopt X, fx, nix, (X/H)q, and px follow- 
ing HKII (listed in Table |2ll. As a metallicity indicator, we adopt 
the oxygen abundance (0 /H), whose sola r value is assumed to be 
12 + log(0/H) = 8.76 i lLoddersll2003h . The abundance of X is 
converted to that of O by assumin g the abundance pattern c alcu- 
lated by the 20 Mq SN model bv Fumeda & Nomotol ( |2002|) (the 
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Figure 1. Grain size distributions adopted in this paper. The grain size dis- 
tributions are scaled with the number density of hydrogen nuclei and the 
metallicity. The solid and dashed lines show the size distributions without 
and with shattering, respectively. The thick and thin lines represent silicate 
and carbonaceous species, respectively. 



Table 2. Adopted quantities. 



Species 


X 


fx 


mx (amu) 


{X/H)0 


px (g/cm^) 


Silicate 


Si 


0.166 


28.1 


4.07 X 10-5 


3.3 


Carbonaceous 


C 


1 


12 


2.88 X 10"'' 


2.26 



ejected masses of O, C, and Si are 1.56, 0.257, and 0.257 Mq, 

respectively). 

In the SN model adopted iUmeda & Nomotol I2OO2I : 
iNozawa et al.l l2007l) . 7^ = 8.87 x 10"^ Vsi = 9.22 x lO"*, 
Vc = 1.71 X 10-^ /in.si = 0.351 and /in,c = 0.135. Besides, 
for the lifetime of a molec ular cloud, we adopt Tci = 10^ yr (e.g. 
iFukui & Kawamurj[2O10h . but we also mention a possibility of a 
longer lifetime in Section l42l 



3 RESULTS 

3.1 Without Shattering 

We examine the relation between dust-to-gas ratio and metallicity 
for the grain size distribution without shattering. The result is ex- 
hibited in Fig.|2}i. At low metallicities [12 + log(0/H) < 8.5], the 
dust is simply supplied from SNe, which leads to an approximately 
linear increase of the dust-to-gas ratio in terms of the metallicity. 
When Vx reaches ~ /in,xJ'x//?SN (~ 3 x 10^ and 2 x 10"^ 
for silicate and carbonaceous dust, respectively), the dust destruc- 
tion becomes comparable to the dust supply from SNe (equation 
Is}. Thus, the increase of the dust-to-gas ratio saturates around 
12 + log(0/H) ~ 8.5. Around 12 + log(0/H) ~ 9 (~ 2 Zq) 
the dust-to-gas ratio rapidly increases because of the nonlinearity 
of grain growth (dDx/dZx oc VxZx', that is, grain growth occurs 
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Table 1. Moments. 







Silicate 






Carbonaceous 






(a)o {^lra) 


(a?)Q (/im^) 




(a)o (Aim) 


(a?)o (miii^) 


(a3)o (^m3) 


Nozawa et al. (2007) 
Hirashita et al. (2010) 


2.32 X 10-2 
7.84 X 10-* 


1.95 X 10-3 
3.23 X lO"*^ 


4.83 X 10""' 
4.35 X lO-'^ 


8.59 X 10-3 
5.42 X 10-* 


3.45 X lO-"* 
5.56 X 10-^ 


6.35 X 10-5 
1.60 X 10-** 



as a result of collisions between dust grains and metals) in molec- 
ular clouds. As shown in HKl 1, the grain growth becomes promi- 
nent at a certain critical metallicity depending on dust grain size 
distribution (see Section |4Tt . Under the size distribution without 
shattering, dust growth becomes prominent only after the metallic- 
ity reaches a super-solar value. This is because of inefficient grain 
growth caused by the large surface-to-volume ratio of dust grains. 

In addition, SNe produce a large amount of silicate than car- 
bonaceous dust. After grain growth in molecular clouds, the silicate 
species is still dominant over the carbonaceous species. 

3.2 With Shattering 

If we adopt the grain size distribution with shattering, we obtain the 
result exhibited in Fig. |2j). Compared with the case in the previous 
subsection, the rapid increase of dust-to-gas ratio occurs at a lower 
metallicity, 12 + log(0/H) ~ 8.3 (~ 1/3 Zq), because shattering 
remarkably raises the surface-to- volume ratio of the dust grains. 

3.3 Observational Data 



iMichalowski et al.U20IOl) have estimated the dust and gas mass by 
using the millimeter continuum and CO lines, respectively, for nine 
> 5 quasars. Their samples show high dust-to-gas ratios from 
1.65 X 10^"^ to 4.30 X lO^'^, which could be regarded as upper 
limits because they did not take neutral hydrogen into account. If 
we alternatively utilize their dynamical mass as an upper limit for 
the gas mass, we then obtain V > 8.0 x 10-* (assuming an in- 
clination angle of 4 0° ) to be a lower limit of the dust-to-gas ratio. 
ljuarez et al. claim that high-z {z > 5) quasars are likely to 

have super-solar metallicities (< 7 Zq) based on observations of 
broad line regions. Narrow line regions (NLRs) may trace the metal 
content of the entire QSO better, and Matsuoka et al. (2009) indi- 
cate that the NLRs in quasars at 1 < 2; < 4 have solar metallicities. 
Considering that the bright QSOs at z > 5 are highly biased to 
the most luminous and metal-enriched systems, it is reasonable to 
assume that the metallicities of 2 > 5 QSOs is > Zq. Comparing 
our results with the observational lower limit of dust-to-gas ratio 
(Fig.|2j, we find that at least 12 + log(0/H) = 9.2 (Z ~ 3 Zq) is 
required to explain the lower-limit of dust-to-gas ratio if we adopt 
the grain size distribution without shattering. Thus, a super-solar 
metallicity is strongly required to explain the rich dust content in 
high-redshift quasars. With shattering, however, the dust-to-gas ra- 
tio reaches the observational lower limit at a more moderate metal- 
licity (12 + log(0/H) = 8.6; Z ~ 0.7 Zq) because of the grain 
growth at a lower metallicity. 



ity is larger than a certain critical metallicity. According to HKl I, 
the critical metallicity Zcr is the metallicity such that I3x{Zct) = 
e/?SN ~ 0.965. For the dust grain size distribution without shatter- 
ing, the critical metallicity of the silicate species and carbonaceous 
species are 12 + log(0/H) = 9.08 and 12 + log(0/H) = 9.35 
(~ 2-4 Zq), respectively. These metallicities actually correspond 
to the metallicity levels where a strong rise of dust-to-gas ratio 
occurs by grain growth (Fig. |2). In contrast, for the dust grain 
size distribution with shattering, the critical metallicity becomes 
12 + log(0/H) = 8.26 and 12 + log(0/H) = 8.25 (~ 1/3 Zq 
for silicate and carbonaceous species. The variation in the critical 
metallicity between the grain size distributions with and without 
shattering explains the difference of the metallicity level at which 
the dust-to-gas ratio is raised by grain growth (see Fig. |2). 



4.2 Dependence on other parameters 

According to equation (|4j, the hydrogen number density in molec- 
ular clouds, nn, is also an important quantity that determines the 
time-scale of grain growth, r. We have presumed that wh = 
10^ cm-3. However, in some environments such as circumnu- 
clear starburst regions in high- 2; quasars, the physical conditions 
of molecular clouds can be quite different. Some recent observa- 
tions imply that the density in mole cular clouds in high-2 quasars is 
higher than that in nearby galaxies ( Riechers et al ]|2007|J Gap et al.l 



l2007l :lKi essen. Spaans. & JaDPsenll2007h . In order to investigate a 
possibility of a larger tih, we adopt riH = 10* cm-^. In Fig.[2};, we 
show the result for the dust grain size distribution without shatter- 
ing. We observe that higher gas density enhances grain growth and 
that the rapid increase of dust-to-gas ratio occurs at a lower metal- 
licity. Therefore, the grain size distribution without shattering can 
be reconciled with the observed dust-to-gas ratio in high-z quasars 
if we assume such a high density as nu > 10* cm-^. 

Equation (3) offers another possibility of enhancing grain 
growth. The dust mass increment /3x is a function of Td/r. As a re- 
sult, a longer lifetime of molecular clouds hosting the grain growth 
Tci has the same effect as a shorter typical time-scale of grain 
growth T, which means that Td = 10** yr instead of Td = lO'^ yr 
has the same eff ect as applying n n ~ 10* cm-^ instead of 
TiH ~ 10 cm' " MKodaetal] j2009l) have indicated that the exis- 
tence of sheared structure within the spiral arms in the Milky Way 
supports a lifetime of molecular clouds comparable to galactic rota- 
tional time-scales (~ 10** yr). Moreover, the gas temperature Tgas 
and the sticking probability Sx also enter the time-scale (see equa- 
tions [6] and |7j; however, since Tgas and Sx cannot be much larger 
than the adopted values, it is difficult to make r shorter by changing 
Tgas or Sx- 



4 DISCUSSION 

4.1 Critical metallicity for grain growth 

llnoud ( I2OIII) and lAsano et al.l ilOlA have shown that the grain 
growth by accretion dominates the grain abundance if the metallic- 



5 CONCLUSION 

The importance of grain size distribution on dust enrichment in 
high-redshift (z > 5) quasars has been investigated in this work. 
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12 + log(0/H) 



We have found that, if the grain size distribution of SN-dust is con- 
sidered, very high metallicity (> 3Zq) is required to explain the 
dust-to-gas ratio of those quasars under a similar quantities to the 
local molecular clouds (tih = 10^ cm~^ and Td = 10*^ yr). How- 
ever, if small grains are produced by shattering, grain growth is ac- 
tivated and the dust-to-gas ratio in high- 2; quasars can be explained 
by more moderate metallicities(> 0.7 Zq). Thus, we conclude that 
the dust grain size distribution has a dramatic impact on the evolu- 
tion of dust content in high-z quasars. 
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Figure 2. Relation between dust-to-gas ratio and metallicity. The solid lines 
shows the total dust-to-gas ratio. The dotted and dashed lines present the 
contributions from silicate and carbonaceous dust, respectively. The hor- 
izontal dot-dashed line show the lower limit for the dust-to-gas ratio in 
the high-2 quasar sample in Michalowski et al. (2010). Panels (a) and (b) 
present the results for the grain size distribution without and with shatter- 
ing, respectively, along with triangles marking the critical metallicities for 
grain growth. Panel (c) shows the result for a high density (njj = 10^ 
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